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ABSTRACT 

The cosmological evolution of Active Galactic Nuclei (AGN) is important for understanding 
the mechanism of accretion onto supermassive black holes, and the related evolution of the 
host galaxy. In this work, we include objects with very low Eddington ratio (10~ 3 — 10~ 2 ) 
in an evolution scenario, and compare the results with the observed local distribution of black 
holes. We test several possibilities for the AGN population, considering obscuration and de- 
pendence with luminosity, and investigate the role of the Eddington ratio A and radiative 
accretion efficiency e on the shape of the evolved mass function. We find that three distinct 
populations of AGN can evolve with a wider parameter range than is usually considered, and 
still be consistent with the local mass function. In general, the black holes in our solutions 
are spinning rapidly. Taking fixed values for e and A neither provides a full knowledge of the 
evolution mechanism nor is consistent with the existence of low Eddington ratio objects. 
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Understanding how Active Galactic Nuclei (AG N) work has b een 
a major goal since the discovery of quasars JSchmidt|[l963h . It 
is now believed that mass accretion onto a supe rmassive black 
hole (SMBH) is the mechanism that powers AGN JSalpete il l 1964 
'. lLvnden-Belllll969l : lReeslll984l) , and this hypothesis is supported 
by recent evidenc e for massive black holes in the centre of most 
nearb y galaxies iKormendv & R ichstone 19931; iRichstone et al.l 
1998). The mass of t he black hole (BH) is related to the propertie s 
of the host galaxy jMagorrian et al.| [l998: Marc oni & Huntll2003h . 
implying a close relation between the evolution of galaxies and the 
BH in their centres . 

Soltan (1982), developed a method to obtain the total mass 
density of SMBH, pbh, from the emitted light and number 
counts of AGN, providing a way to trace the accretion history 
of these BHs, from the AGN luminosity function as a function 
of redshift. Following this approach, there have been attempts 
to compute the evolution of the BH mass function from a start- 
ing redshift, and compare it with the BH distri bution in nearby 
galax i es, kno wn as the local distribution JSmall & Blandfordl 



19921; ISalucci et al.l [l999fc iMarconi et all 
Shankar. Weinberg. & Miralda-Es cude 2007). Progress has been 
made by the use of semi-analytical models and increasingly bet- 
ter observational results, although it is still hard to constrain the 
important factors affecting this calculation: the luminosity function 
at each redshift, the bolometric correction, the obscuration and the 
free parameters related with the accretion process. 

There is evidence from the X-ray background that obscura- 
tion plays an important role in AGN evolution. The background 
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is due to the in tegrated emission from A GN with high intrinsic 
absorption (e.g. | Bran dt & Hasinger 2005), and by modelling it, 
iFabian & Iwasawal J 19991) found that most of the accretion onto 
AGN is obscured. One of the main goals of AGN studies is to 
get unbiased samples of objects, and for this, X-ray surveys reveal 
more AGN than found in other wavelength bands. The observations 
in the hard X-ray band (2-10 keV) are the least biased against ob- 
scuration, and th erefore more suitable to observe and study AGN. 
lUeda et alj2003l (hereafter U03) show that the hard X-ray luminos- 
ity function is best-fitted by a luminosity dependent density evolu- 
tion (LDDE), and that the number density of low luminosity AGN 
peaks at lower redshift than the high luminosity ones. To be able to 
convert the X-ray to total energy output, and obtain a knowledge of 
the multi-wavelength properties, a bolometri c correction is needed , 
this is a function of source lumino sity (e.g.. [Marconi et alj|2004l) . 
or more likely, the Eddington ratio dVasudevan & Fabia n 2007). 

The accretion mechanism is another important factor. There 
are two important free parameters to describe it; one, the radiative 
accretion efficiency e, representing the fraction of accreted mass 
radiated, establishes the relation between the bolometric luminos- 
ity and the growth rate of the BH, the other, the Eddington ratio A, 
relates the AGN bolometric luminosity with the Eddington lumi- 
nosity: L = AZ/e, with Z/e = 4ncGMrriH /<tt . In previous work, 
the Eddington ratio has been taken to be either a cons tant, changing 
as a f unction of redshift, or a function of mass (see lShankar et alj 
2007 and references therein). 

Another factor is the hydrogen column density (A*h) dis- 
tribution, often called the A^h function. Several studies tried to 
parametrise it, U03 showed that the absorption is luminosity depen- 
dent, with A^h a decreasing function of the intrinsic luminosity, and 
found no significant dependence on redshift, although this question 
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is stil l in debate (see for example lTozzi et al j2006l : lTreister & Urrvl 
l2006h . 

To constrain the parameters involved in the growth of AGN, 
Marco ni et al, I J2004h compare the mass function from AGN relics 
with the local mass function, using a single AGN population. They 
find an acceptable match for 0.1 < A < 1.7, and 0.04 < e < 0.16, 
with the values e = 0.08 and A = 0.5 giving the minimum 
deviation between the two mass functions. However, recent re- 
sults showed that the Eddington ratio can be si gnificantly lower , 
with many objects having A ~ 10~ 2 - 10" 3 , teabic et al.ll2007l ; 



Fabia n. Vasudevan. & G andhi 2008), at least one order of magni- 
tude less than the Marco ni et al 1 120041) lowest acceptable value. 
This indicates that there is a wide spread in the range of A, and that 
we have to consider low values for this parameter when computing 
the AGN mass function. 

In this paper, we analyse how the free parameters A and e in- 
fluence the shape of the evolved AGN mass function, starting by 
considering one population of AGN, which is then split in two, ac- 
cording to obscuration. This leads to an implausible excess of high 
mass objects so we add a luminosity dependence for AGN with low 
A, and consider the evolution of three different populations. The 
real situation is undoubtedly more complex than our simple mod- 
els, but we expect to have captured the essence of AGN evolution 
in the (e, A) plane. 

We adopt standard cosmological parameters of Ho = 70 km 
s" 1 Mpc" 1 , Cl m = 0.3 and fi A = 0.7. 



2 THE BLACK HOLE MASS FUNCTION 

To obtain a black hole mass function for z — 0, and compare it 
w ith the local observed mass function, we follow the procedur e 
of lMarconi etai] J2004I) (hereafter M04) and lShankar et all fe007l) . 
Starting with an AGN luminosity function (LF) at certain redshift 
Zi, and considering that all BHs were active, we follow their evolu- 
tion with cosmic time, assuming that their growth and energy out- 
put are due to mass accretion only. 

Using Soltan's argument (1982), we can establish a relation 
between the total energy output, given by the integrated luminosity 
function, and the density of BH: 



PBH 



i-6 /•*> dt r L 
c 2 inioi Tz z J Ll 



d<t>{L,z) 
d log L 
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Here, L is the bolometric luminosity, e the radiative accretion effi- 
ciency and: 

£ = - [(1 + z)H i) ^(l + zY i Q ln + n x ] _1 , 

for the cosmological model used in this work. We are interested 
in the AGN distribution of masses, i.e., what is the number of ob- 
jects with a certain mass M, per unit of comoving volume, and 
how it varies with redshift. Fo r this, we use the continuity equation 
for the non-merging case , (see[Cavaliere, Morrison. & Woo dll97ll : 
ISmall & Blandforj|l992h . Considering that the BH growth is only 
due to mass accretion, we assume that the total number of BHs is 
constant with time, and obtain: 



d d(f>(M,t) 
dt dM 
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(1) 



In the case of merging, the right side of the continuity equation 
would not be zero, but equal to a source function, implying that 
the total number of BHs is not conserved. In Eq.Q] <f>(M, t) is the 
number of BH with mass M per unit comoving volume for a certain 



cosmic time t and (M(M, t)) is the mean accretion rate for BH of 
mass M (active and inactive): 



(M(M,t)) = 5(M,t)M(M,t). 



(2) 



Here S(M,t) is the fraction of active (accreting) BH as a function 
of M and t, called the duty cycle. M(M, t) is the accretion rate (or 
growth rate) of a BH with mass M at time t. 

The growth rate is a result o f the mass that falls t o the BH but 
is not converted into energy (see iMarconi et al! I l2004l) . If we have 
a total mass falling (A/ to t), part of it is converted into energy L = 
eMtotC 2 , and the rest will accrete onto the BH: M = (1 — e)Mtot- 
Consequently, 



L = 



(1 



M(M, t). 



The energy output at a certain t is due to active BHs, so 

d ^ t l d lo S L = S(M,t) dm ^dM. 

dM 



dlogL " ° -v— dM —• (3) 

Using Eq.[2]and Eq.[3]to eliminate 8(M,t), and replacing in the 
continuity equation Q3 
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= 0. 



— , and as- 



suming constant Eddington ratio (A) and e, one can solve the con- 
tinuity equation to obtain the mass function at a redshift Zf. 

d4>(M,Zf) _ (1 - e)A 2 c 2 



dM 



et| In 10 



f z ' d 
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L t dL 
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dz 
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(4) 



This equation tells us that the mass function at a certain redshift, 
can be obtained with a Soltan-type argument, by tracking the en- 
ergy that was produced due to mass accretion from zi to Zf. We 
assume that the initial mass function (at redshift 3), does not con- 
tribute significantly to the final mass function, since the major black 
hole growth occurs for z < 3 (M04). 

To be able to compare this evolved mass function with the lo- 
cal one, the previous equation has to be evaluated for Zf = 0. We 
have to assume a certain luminosity function, a bolometric correc- 
tion, set the initial conditions and give values to e and A. For the LF, 
we consider the luminosity-dependent density evolution (LDDE) of 
U03, which describes 9 !?l Lyi T ' z ' > , the number of AGN with intrinsic 
luminosity Lx per unit of comoving volume per logLx, as func- 
tion of the luminosity and the redshift: 



d<j)(L 



x,z) 



d logLx 
where e(z. 

e(z,Lx) = 



Lx 
L t 



+ 



Lx 
Lt 



Lx) is the evolution term, 

(1 + 2) P1 



e(z, Lx) 



z < z c {Lx), 



Zc(Lx) = 



[1 + z c (L x )] pl - p2 (1 + z)^ z > z c (Lx), 
and the redshift cut off, 

Lx > L a , 

i^y ^<L a . 

The parameters have the values: A = 5.04 x 10 -6 Mpc -3 ; L a = 
10 44 6 erg s _1 ; L* = 10 43 94 erg s" 1 ; 7l = 0.86; j 2 = 2.23; pi = 
4.23; p2 = -1.5; z* = 1.9; a = 0.335. This luminosity function has a 
discontinuity in its derivative, which leads to unphysical solutions 
for the mass function. To avoid that, we rewrite the luminosity func- 
tion to make it continuous in its derivatives (A. Marconi, private 
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communication) : 



d<j)(Lx, t) 
d log L x 



^ dd\ • L. : . 



d log L> 



+ 



( d(j>(L x < L a ,t) 
\ d log L x 



with /3 = 8. 

For the bolometric luminosity L, we use the Edding- 
ton ratio dependent bolom etric correction k(A) = L/Lx 
dVasudevan & Fabianl l2007h . with k(X) = 19.3 for A < 0.1, 
ft(A) = 54.5 for A > 0.3, and an intermediate value of re(A) = 
36.9 for the other Eddington ratios. To account for Compton- 
thick AGN, we multiply our mass function by a c onstant factor 
of 1.6, following the result by Risaliti et al] dl999h . In this work 
we consider Z\ = 3, and that all BHs were active at that time 
(5(M, Zi ) = 1). Using 



dlogL 



d(j)(Lx, z) 



dlog L x , 



d(j){L,z) 

d log L ° d log Lx 

and since k(A) is constant for an assumed value of A, we finally 
obtain 

OL 
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For a simple model of one population of AGN, we apply the method 
described before, noting that each set of parameters A and e give a 
different mass function for the same luminosity. To add the effect 
of obscuration, and understand the role of the absorbed AGN in the 
shape of the mass function, we consider the AGN population to be 
separated in two groups with different absorption column densities 



JVh: unabsorbed (10 
(10 22 cm~" 



cm -2 < JVh < 10 . 
< Nn < 10 24 cm- 2 ). We define 



ATh(Lx) 



N-2 



/(Lx;Af H )dlogAr H , 



), and absorbed 



(6) 



where f(Lx',Nn) is the probability distribution function from 
U03, normalised to unity in the Compton-thin region 20.0 < 
log/V H < 24.0. We find A/"n bs (Lx) integrating Eq.|6] between 
the values of column density assumed for absorbed AGN, and 
A/n n (I/x) integrating in a similar way for the unabsorbed ones. 
That is, for each value of Lx, we split our population in two groups, 
with constant number fractions given by the J\fn(Lx) function. In 
summary we have: 



d<j>(M, 0) 1 
dM 

d4(M, 0) ' 
dM 



■ M^iLx) x 



d<t>(M, 0) 
dM '■ 

d<f){M, 0) 
dM ' 



(7) 



(8) 



3 CONSTRAINTS ON THE FREE PARAMETERS 

To constrain the free parameters A and e, we compare our evolved 
AGN mass function with the local mass function. As a reference, 
we use the local mass function in M04 and the respective la er- 
ror, in the range (10 6 - 10 10 )M Q which was obtained using both 
the Mbh — a* (velocity dispersion), and the A/bh — Lbui (bulge 
luminosity) relations. The method to obtain the mass function is 
as follows: for each mass value, we compute the corresponding 
Lx by first using L — ALe to obtain the bolometric luminosity 
from the mass, and then k(A) to convert from bolometric to intrin- 
sic luminosity. We then replace the U03 luminosity function and 
the bolometric correction in Eq. [5] differentiate in L and integrate 



— 



10 



log M B11 [M G 



Figure 1. Comparison between the evolved nu mber of AGN per un it vol- 
ume and the local mass function (in black) by Marconi et al. ( 20q3). The 
grey area is the la error in the M04 local mass function. The green curve 
represents the unabsorbed AGN, and the blue curve the absorbed AGN. 
The red curve is the total population (absorbed plus unabsorbed). The val- 
ues used were A„hs = 10 ~ 3 , e a b s = 0.085, A un = 0.56, e un = 0.21 
with the lUedaet aU 120031) Nn distribution. For a low value of A, a two 
population model predicts a higher number of high mass black holes than 
that observed. 



from z = to z = 3 to obtain ■ For the case of two or 

more populations we use Eq. [6] to compute the A/h function for 
each value of Lx, and replace it in Eq.|7]for the absorbed and Eq. 
[8] for the unabsorbed, allowing each population to have a different 
Eddington ratio and accretion efficiency. 

For one population, with only two free parameters, we find no 
value of e or A for which the evolved mass function is consistent 
with the local one within a la error. Considering a black hole mass 
range of (10 r — 10 9 )A/q, we obtain agreement but only for high 
Eddington ratios (A > 0.1). To search for alternative scenarios and 
be able to have lower A we study the effect of obscuration, by con- 
sidering the A^h function by U03. 



3.1 Two populations 

We start by assuming two populations, one of obscured and one 
of unobscured AGN. The agreement between the local and the 
evolved number density of AGN can be obtained by changing si- 
multaneously the free parameters for the absorbed and unabsorbed 
objects. To constrain the allowed ranges for the accretion effi- 
ciency and the Eddington ratio, we search parameter space, find 
the evolved number density of AGN for every combination of the 
four parameters (A a b s , e a bs, A U n, tun), and compare it with the lo- 
cal number density from M04. We search the parameter range: 
1(T 3 < A < 1, 10~ 2 < e < 0.42. The upper value of e is limited 
by the maximum accretion onto a Kerr black hole, if we consider 
a case of a radiatively efficient system, the radiative accretion ef- 
ficiency e is equal to the mass to energy efficiency (rf). In a more 
general case, e < rj. 

With two populations, we do not find parameters that re- 
sult in a mass function within la of the local mass function. We 
only find good agreement if we limit the mass range to Mbh = 
(10 7 — 10 9 )Mq, but for low values of A, it predicts a higher density 
of high mass black holes (> 10 9 Mq) than that observed, which is 
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not consistent. In Fig. 1 we chose a set of values to illustrate our 
problem, the total mass function in red is compared with the lo- 
cal mass function from M04. The blue and green lines represent 
the absorbed and unabsorbed populations respectively, with fixed 
parameters. For low values of A, although it is possible to obtain 
good agreement in a limited mass range, the sum of the absorbed 
and unabsorbed populations (red line), predicts a higher number 
of high mass black holes than observed. Thi s is a common prob- 
lem, b oth for M04 luminosity dependent and lVasudevan & Fabianl 
d2007h Eddington ratio dependent bolometric corrections. When A 
is lowered (A < 1CP 2 ), the mass function shifts down and right, 
which means that it raises the number density of high mass black 
holes. Nevertheless, low A objects are observed in the Universe, 
and we require that they agree with the local mass function. 



3.2 Three populations 

To obtain agreement in the high mass black hole range, the num- 
ber of AGN with low Eddington ratio has to drop with luminosity. 
By hypothesis, and in agreement with the TVh function in U03, we 
assume that the low A objects are included in the obscured pop- 
ulation. We therefore divide our absorbed population in two: the 
absorbed objects with lower Eddington ratio (10~ 3 < A < 10 -2 ), 
and objects with higher Eddington ratio (10~ 2 < A < 1). This 
way, our AGN population is divided into three groups: obscured 
objects with lower A and numbers dropping with increasing lu- 
minosity, obscured objects with higher A, and unobscured objects 
(1(T 2 < A < 1). 

For the luminosity dependence, we define a simple decay for 
the fraction /(Lx), the number of objects with low A divided by 
the total number of absorbed objects for a certain luminosity: 

/(Lx) = 51og(L x )+6, 

where S is the free parameter, and b is determined by assuming that 
for Lx < 10 41 erg s" 1 all the absorbed AGN have low A (/ = 1). 
Therefore, the mass functions for the three populations are: 

- Absorbed with low A: 

#(M,0) ' abS WA/ *.,r N .. #(M,0) 

— — =/(- L x)x7V H (Lx) x 

- Absorbed: 

d<t>{M, 0) abs 
dM 

- Unabsorbed 
#(M,0) un 



dM 



= (1 - f(Lx)) x A^ bs (L x ) x 



d4(M, 0) 
dM 



= MS 1 (Lx)x 



dcj>(M, 0) 



(9) 



(10) 



(11) 



dM K ' dM 

In terms of free parameters, we now have 
(A a bs, e a bs, A un , e m , Ai abSj ei abB , S). Following this method, 
we find several sets of parameters for which there is good agree- 
ment with the local mass function. In Fig. 2 we have an example 
of a set of parameters for which there is good agreement, the grey 
area is the lcr error of the M04 local mass function (in black) 
and the red line is our total evolved mass function. The green 
line represents the unabsorbed population, and the blue lines the 
absorbed: light blue for the low A and dark blue for the others. 
As we can see, a decay in / with luminosity allows the low A 
population to agree in the high black hole mass range, and we 
obtain very good agreement between our evolved mass function 
and the local one. In terms of the effect of a parameter change, 
lowering e, makes the mass growth of the BHs more efficient, 
and the curves in Fig. 2 go up. A shifts the curves to the high 
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Figure 2. Comparison between the evolved nu mber of AGN per u nit vol- 
ume and the local mass function (in black) by Marconi et al. ( 20oJ. The 
grey area is the lcr error in the M04 local mass function. The green curve 
represents the unabsorbed AGN, and the blue curves the absorbed AGN: 
light blue for the absorbed objects with low A and dark blue for the others. 
The red curve is the total population (absorbed, absorbed low A and unab- 
sorbed). The values used were Ai a j, s = 0.007, ei a b s = 0.21, A a i, s = 0.5, 
e a .h R = - 068, A un = 0.079, e un = 0.21 and S = -0.3 with the 
lUeda et al l 120031) Nn distribution. 



mass or low mass end. An increase in A corresponds to a shift 
to the left and up, a decrease makes the curves go right and 
down. The plot in Fig. 3 shows the result of our search in the 
parameter range described for the two populations and also with 
1(T 3 < Ai a bs < 10~ 2 , 1(T 2 < eiabs < 0.42 and -1< S < 0. 
Each population is described by a set of three points in the plot: 
light blue (light grey), dark blue (black) and green (medium grey). 
The sum of these three populations gives the total one, which is 
evolved and compared with the local mass function. We plot all 
the values of A and e for which the evolved mass function is within 
a lcr error from the local one. Note that the number of points are 
not the same for all populations because an unabsorbed population 
can have, for example, more than one absorbed or low lambda 
absorbed population for which the mass function agrees well with 
the local one. 

In terms of the general features, we see that we obtain a wide 
range of possible parameters. Each population is described by a 
point in Fig. 3, which means that the total mass function (red line 
on the previous plots), is described by a definite set of three points, 
light blue (light grey), dark blue (black) and green (medium grey), 
and not a random combination. The absorbed and unabsorbed, blue 
(light grey and black) and green (black), populations, have val- 
ues of A that are in general higher that 0.1, and the absorbed with 
low A, light blue (light grey), spread in the whole parameter range 
(10~ 3 - 10" 2 ). We notice a void, due to the fact that for the ab- 
sorbed with high A and unabsorbed, an Eddington ratio > 0.1 
causes the evolved mass function to disagree from the local mass 
function in the high mass range. Allowing the Eddington ratio for 



the low A absorbed objects to go higher than 10~ would fill in the 
void. Values for the slope S, are found between —0.75 and —0.25, 
which means that the values for the luminosity Lx where / reaches 
zero range from 10 42 ' 3 erg s _1 to 10 45 erg s _1 . These values corre- 
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spond to a steeper de ca y than for example the absorbed fractions 
from lHasingerl J2008h or lTreister et al.l d2008h . which combine ab- 
sorbed objects of all A. The values for the accretion efficiency are 
above 0.035, and in this scenario, lower for absorbed objects than 
unabsorbed objects in the same range of A. We do not find evidence 
for advection dominated accretion flows (ADAF) behaviour in our 
low efficiency objects. As we can see from Fig 3, almost all possi- 
ble parameters are above e = 0.057, which is the efficiency corre- 
sponding to a Schwarzschild (non-rotating) black hole, indicating 
that the black holes must be spinning. In fact, most of the black 
holes in our solutions are rapidly spinning. A population of low A 
objects with low efficiency, light blue (light grey) points, or low 
efficiency absorbed objects, dark blue (black) points, can only rep- 
resent the local mass function if combined with unabsorbed, green 
(medium grey), objects with high efficiencies (e > 0.1). 

We have assumed that Compton-thick objects contribute to the 
mass function equally in all the mass range, and set a correction fac- 
tor of 1.6 to account for them. A change in this value would not af- 
fect qualitatively the shape of the relation between e and A showed 
in Fig. 3. To evaluate the errors, we used two approaches: a higher 
multip lying correction factor and a similar method to iGilli et al.l 
(2007). In terms of the mass functions, a higher multiplying fac- 
tor shifts all the curves up, making the lower values of e in Fig. 3 
be inconsistent with the local mass function. In summary, a higher 
fraction of Compton-thick objects would imply that the local mass 
function can be reproduced using higher a ccretio n efficiencies than 
the ones from Fig. 3. For the IGilli et al.l J2007I) Ah function, we 
obtain results that do not differ significantly from the ones with the 
U03 Ah function, the main difference is that the parameter values 
for the absorbed objects (high A), are found in the same range as 
the unabsorbed, that is, higher than the ones from U03. 



4 SUMMARY AND CONCLUSIONS 

We have tried to reconcile the existence of low Eddington ratio 
AGN, with the local mass function. Our approach has been to try 
the simplest solution and then add features to obtain good agree- 
ment. In order to include a population of low A objects and a 
large spread in possible A values, we find that at least three pop- 
ulations are required, with different absorption properties depend- 
ing on intrinsic luminosity for the lowest A objects. AGN evolu- 
tion is not a straightforward problem, and may in future require 
more complex behaviour to be considered. We find that most of 
the black holes in our solutions have high accretion efficiencies 
(e > 0.1), and therefore must be spinning rap idly, which has impli- 
cations for the merger and accretion histories dVolonteri et al.l20o3 ; 
iKing & Prinriell2006t) . 
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